Oscillons After Inflation 



O 

> 
O 

(N 

o 

u 

43 

Oh 
i 

o 

CO 

> 

in 

m 
m 
m 

O 



> 

X 



Mustafa A. Amin 1 , Richard Easther 2 , Hal Finkel 2 , Raphael Flauger 2 , and Mark P. Hertzberg 3 
1 Department of Physics, Massachusetts Institute of Technology, Cambridge MA 02139 

2 Department of Physics, Yale University, New Haven CT 06520 and 
3 SITP and KIP AC, Department of Physics, Stanford University, Stanford CA 94305 

Oscillons are massive, long-lived, localized excitations of a scalar field. We show that in a class of 
well-motivated single-field models, inflation is followed by self-resonance, leading to copious oscillon 
generation and a lengthy period of oscillon domination. These models are characterized by an 
inflaton potential which has a quadratic minimum and is shallower than quadratic away from the 
minimum. This set includes both string monodromy models and a class of supergravity inspired 
scenarios, and is in good agreement with the current central values of the concordance cosmology 
parameters. We assume that the inflaton is weakly coupled to other fields, so as not to quickly 
drain energy from the oscillons or prevent them from forming. An oscillon-dominated universe has 
a greatly enhanced primordial power spectrum on very small scales relative to that seen with a 
quadratic potential, possibly leading to novel gravitational effects in the early universe. 



Simple, single-field models of inflationary cosmology 
are often associated with energy scales far beyond the 
reach of present day accelerators and the properties of 
the post-inflationary universe are largely unknown. One 
constraint on this phase is that energy must be extracted 
from the oscillating inflaton condensate, ensuring that 
the universe becomes radiation dominated, setting the 
scene for the hot big bang, and the production of the 
cosmological neutrino background and nucleosynthesis. 
A widely-studied candidate for this process is paramet- 
ric resonance [THS]. In many cases, the potential can be 
se//-resonant, where resonance generates quanta of the 
inflaton field itself more efficiently than particles cou- 
pled to the inflaton. We show that in a class of well- 
motivated self-resonant models, the universe may become 
dominated by oscillons: massive, localized, metastable 
configurations of a scalar field [4l413j. 

In this letter, we study a single inflaton, 0, with a 
canonical kinetic term and potential, V(<p), minimally 
coupled to Einstein gravity. Oscillons can form if 
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where U(0) = and U(4>) < for some range of <j> (see, 
e.g. [TO])- Consider potentials with V(<p) ~ (j) 2a during 
inflation and a < 1. These are generated by a number 
of string and supergravity scenarios |14H20j . and yield 
U{4>) < at large <j>. We require that V{<j>) has a stable 
minimum which we chose to be at the origin, so it is 
natural to expect that V(cf>) ~ 4> 2 for small 4>. Finally, 
by continuity, there is necessarily some crossover scale, 
(j> « M, between these two regimes. We capture this 
with the following explicit potential 



V(<f>) = 
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M 2 J 
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The precise forms of V(<p) in scenarios with V ~ (j> 2a 
can differ from equation ([2]): however, our results sug- 
gest that while oscillon formation is sensitive to M, it is 



insensitive to the detailed form of the potential. More- 
over, for a = 1/2 we reproduce the axion monodromy 
potential [HI [T5]. We stipulate that the couplings be- 
tween the inflaton and other fields are small enough for 
them to be ignored. 

The tensor-scalar ratio, r, and scale-dependence in the 
scalar perturbations, \n s — 1|, grows with a (see, e.g. 
[2"T]). Quartic inflation (a — 2) is ruled out by cur- 
rent data [2"2"H2"1] and even quadratic inflation (a = 1) 
is somewhat disfavored, relative to models with a < 1 
|25) . Consequently, the above potential is well- motivated, 
both theoretically and phenomcnologically. 1 

The post-inflationary universe is initially smooth, so 
even if a potential supports oscillon solutions, an actual 
oscillon-dominated phase requires a mechanism for gener- 
ating inhomogeneity within the post-inflationary horizon 
(see, e.g. jllj). Equation [2] supports parametric reso- 
nance when a < 1, which lead to the explosive produc- 
tion of 4> quanta, and a highly inhomogeneous universe. 
However, M is large, the V(<j>) effectively quadratic dur- 
ing both the last portion of inflation and subsequent os- 
cillatory phase, suppressing resonance and oscillon pro- 
duction. Conversely, if M is significantly sub-Planckian 
we see resonance and oscillons can form. Note that nar- 
row resonance also occurs when a > 1, but V(<f>) cannot 
support oscillons. 

In what follows, we first summarize the inflationary 
dynamics and describe a Floquet analysis of the resonant 
phase. We show that strong resonance and a subsequent 
oscillon-dominated phase requires < a < 0.9 and M < 
0.05Mpi (M p i = l/v / 87rG r /v), which may be realized in 
the physical scenarios that motivate these models. We 
then discuss the cosmological consequences of an oscillon- 
dominated phase. 



1 Oscillon production in hybrid inflation models is studied in |9]. 
These models have n s > 1, and are disfavored by observations. 
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INFLATIONARY DYNAMICS 



The observed amplitude of the primordial fluctua- 
tions effectively removes one free parameter from the 
potential in equation |2]). Further, we will see that we 
are primarily interested in models where M is substan- 
tially smaller than the Planck mass, so that V(4>) ~ 
m 2 M 2 (4>/M) 2a /2a during inflation. 

Astrophysically interesting perturbations are laid 
down when the remaining number of e-folds before the 
end of inflation, N ~ 55, though in general, N is a func- 
tion of the post-inflationary expansion history |21j . Using 
standard slow-roll approximations, the amplitude of the 
power spectrum of curvature perturbations is 
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(3) 



For a given a and /3 = 
with N = 55 and A 2 R 



M p i/M, we use the above equation 
= 2.4 x 10" 9 [23] to deduce m. 



RESONANCE AND OSCILLONS 

Oscillon production at the end of inflation with 
U(4>) = -A0 4 /4 + 5 2 6 /6m 2 was studied in [TTJ [13] . 
When {X/g) 2 <C 1, oscillons are copiously generated, 
with properties matching analytic predictions |10l lllj . 
However, an inflationary phase where V(</>) is dominated 
by a 4> e term has an unphysical perturbation spectrum. 
If the above U(<p) is viewed as a truncation of eq. ^ 
then (X/g) 2 ~ 1, and we cannot appeal to the results of 
[TUlfTTlUB] for the properties of oscillons. Thus, to study 
oscillon formation in this physically reasonable scenario, 
we rely on numerical simulations. 

We can gain a heuristic understanding of oscillon for- 
mation by looking at the instability diagram for the 
potential in ([2]), as resonance generates large inhomo- 
geneities which then relax to form oscillons [TT] . Ignoring 
expansion and working in the limit where <f> is approxi- 
mately homogeneous, Floquet theory allows us to write 
the individual momentum modes of 6 as 



= P+(t) 



,/j-kt 



P-(t) 



(4) 



where P± (t) are periodic functions and ipk are called 
Floquet exponents. Our first task is to calculate these ex- 
ponents: if the real part of fi^, $t(pk), is nonzero and its 
magnitude is larger than the Hubble parameter, H ~ t , 
at the end of inflation, the mode will grow. Roughly 
speaking, if \$t([j,f.)\/H > 10, we have strong resonance. 
In an expanding universe, has a physical wavenum- 
ber k/a(t) and thus moves through a number of Floquet 
bands as the scale factor, a(t), grows, as shown in Fig- 
ure [TJ For our potential, with (3 = M p \/M, one can 
show that the maximum value of \^t(fik)\/H as the modes 
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FIG. 1: Floquet diagram with a = 1/2, = 100. The stable 
regions are dark red. Within the unstable bands, lighter colors 
correspond to larger real-valued Floquet exponents. White 
lines show k/a(t) for representative modes in an expanding 
universe. 



traverse the Floquet bands is [\$t(p,k)\/H]max w A(a)f3 
where A(a) « (1/2) [(1 - a) - (1/10)(1 - a) 2 ] . 

We studied the nonlinear dynamics of resonance fol- 
lowing inflation driven by eq. ([2| using PSpectRe 
[26] . PSpectRe solves the fully nonlinear three di- 
mensional Klein-Gordon equation in an expanding back- 
ground whose behavior is governed by the usual Fried- 
man equations, sourced by the average density and pres- 
sure. The backreaction of metric perturbations on the 
field is ignored. Our simulations begin at the first in- 
stant = 0, although our results are insensitive to the 
details of this choice. The scalefactor a = 1 at the begin- 
ning of our simulations. We assume a standard spectrum 
of initial vacuum fluctuations, although we checked that 
our results are qualitatively insensitive to the detailed 
form of the initial conditions. We ignore backreaction of 
the metric perturbations on the field evolution - these 
can be shown to be small during resonance. The initial 
boxsize is L = 25 /m with 256 3 points in the (comoving) 
simulation volume. 

A single timeslice of a representative simulation is 
shown in Figure [2] Given that oscillons are large over- 
densities, a necessary condition for oscillon domination 
is that 



p>2(p) 



P dV 



JpdV 



(5) 



the fraction of the total energy density contributed by 
regions where p/(p) > 2, is nontrivial. Oscillons are 
effectively fixed in space, persisting for a Hubble time 
or more, so the overall density in an oscillon-dominated 
universe at different times is strongly correlated. Heuris- 



3 




FIG. 2: Oscillon configuration with a = 1/2 and /3 = 50. The 
top plot shows regions where p/{p) > 4 (transparent) and 12 
(solid), while the lower plot shows p/{p) on a two dimensional 
slice through the simulation. Length units are l/a(t)m, and 
these plots were made when a(t) = 5.46. 

tically, / > 0.3 is sufficient to guarantee that the field 
configuration (and thus the post-inflationary universe as 
a whole) was dominated by oscillons. Figure [3] shows / 
as a function of a and /?, along with the maximal value 
of the resonance parameter / H] max . We see that 

strong resonance, or [|3? {lik)\/H\max ~ 10, is both neces- 
sary and sufficient for prompt, copious oscillon formation. 

In models for which / is non-zero, it remains approxi- 
mately constant for a several Hubble times after the onset 
of oscillon domination, demonstrating that this phase is 
long-lived, relative to prevailing cosmological time scales. 
Unlike the oscillons studied in HH H3] which have a 
stable, radial envelope, $(r), which evolves very slowly 
with time, here the corresponding envelope is a periodic 
function of time, and the oscillon "breathes" in and out. 
The detailed dynamics of these oscillon solutions will be 
discussed in a future publication, but we have simulated 
a single oscillon (ignoring expansion) over a long interval 




FIG. 3: The statistic, /, is shown at a(t) = 7 (a(t) = 1 at 
the beginning of the simulation) as a function of a and j3 = 
Mpi/M. Contours show maximal value of the [|9?(/ifc)|/i?] max . 
The thick black contour denotes [\$t(p,k) \ / H] max = 7 whereas 
the thin white ones correspond to [\$l(pk)\/ H] ma x = 1,3. 

for representative values of a and j3 after imposing strict 
radial symmetry, reducing the problem to a 1+1 PDE. 
Even though oscillons are not protected by a conserved 
charge and radiate energy J^S-fHj , these simulations sug- 
gest that they live long enough for the universe to grow 
by a factor of 100 or more, and we expect this to be true 
even if the assumption of radial symmetry is dropped. 
Also, the quantum radiation will be small in the regime 
where the self couplings, such as A ~ m 2 /M 2 , are small 

CONSEQUENCES AND DISCUSSION 

We have demonstrated that for a large class of mod- 
els in excellent agreement with the current concordance 
cosmology inflation is naturally followed by an oscillon- 
dominated phase, provided that the couplings to other 
fields are small. These oscillons are generated by para- 
metric resonance, which occurs if the inflationary poten- 
tial turns over from the slow-roll regime to a quadratic 
regime at a scale M <C M p \. 

The inflationary models here are self-resonant, so os- 
cillon production does not require specific couplings to 
other fields. It is likely that any significant couplings 
between the inflaton and other fields can inhibit the for- 
mation of oscillons, by allowing resonant production of 
quanta of these additional fields. Further, couplings to 
other fields can reduce the stability of oscillons by pro- 
viding an additional channel into which they can radiate 
energy. Lastly, the impact of interactions between oscil- 
lons is largely unexplored (however, see |35 j ) . 

Many resonant models include light fields, leaving the 
universe in an intermediate state between matter and ra- 
diation [3D1 [3T] , but massive self-resonant models lead to 
an oscillon-dominated universe that is effectively mat- 
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ter dominated. Our simulations do not include local 
gravity, but perturbations with sub-horizon wavelengths 
will grow gravitationally during the oscillon-dominated 
phase. The same behavior is seen in non-resonant models 
with an (almost) homogeneous inflaton condensate oscil- 
lating in a pure m 2 </> 2 potential 32J. However, in this 
case the primordial density fluctuations are 0(1O~ 5 ) at 
the scale of the horizon and take a long time to become 
nonlinear. By contrast, fluctuations grow rapidly in a 
self-resonant model, leading to a significant enhancement 
in the primordial power spectrum for high k. The possi- 
bility of gravitational collapse and even primordial black 
hole formation during this phase must be carefully ana- 
lyzed [3ll [34] . Given that the oscillons exist on comov- 
ing scales vastly shorter than those which contribute to 
large-scale structure formation, oscillon formation is un- 
likely to directly modify the primordial power spectrum 
on present-day astrophysical scales. However, for any in- 
flationary model the observed power spectrum is a func- 
tion of the post-inflationary expansion history [3TJ [M] . 
Thus, it will be important to account for the existence 
and duration of any matter-dominated phase, oscillon- 
dominated or otherwise, when computing the detailed 
predictions of the model. 

In summary, we have shown that a significant class of 
realistic inflationary models can naturally lead to copi- 
ous oscillon production following inflation, and that these 
oscillons can - for a time - dominate the overall matter- 
density of the universe. This provides a dramatic exam- 
ple of the potential importance of nonlinear dynamics in 
scalar fields to the properties of the very early universe. 
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